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Based on the connection between the Ising model and a 
correlated percolation model, we calculate the distribution 
function for the fraction (c) of lattice sites in percolating clus- 
ters in subgraphs with n percolating clusters, / n (c), and the 
distribution function for magnetization (m) in subgraphs with 
n percolating clusters, p n (m). We hnd that f n (c) and p n {m) 
have very good hnite-size scaling behavior and they have uni- 
versal finite-size scaling functions for the model on square, 
plane triangular, and honeycomb lattices when aspect ratios 
of these lattices have the proportions l:\/3/2:\/3- The com- 
plex structure of the magnetization distribution function p(m) 
for the system with large aspect ratio could be understood 
from the independent orientations of two or more percolation 
clusters in such system. 

PACS numbers: 05.50.+q, 64.60. Ak, 75.10.-b 



Universality and scaling are two important concepts 
in the modern era of critical phenomena Q and for ana- 
lyzing the simulation or experimental data of finite criti- 
cal systems, one often appeals to finite-size scaling 
where both critical exponents and finite-size scaling func- 
tion play important role. The universality of critical ex- 
ponents was well known for a long time |jj , but the uni- 
versality of finite-size scaling functions received much at- 
tention only in recent years JE] In 1984, Privman and 
Fisher first proposed the idea of universal finite-size scal- 
ing functions (UFSSF's) with nonuniversal metric factors 
f§. In 1995-1996, Hu, Lin and Chen (HLC) ji applied 
a histogram Monte Carlo simulation method |ll]] to cal- 
culate the existence probability [jfl) (also called crossing 
probability |12|) E p , the percolation probability P, the 
probability for the appearance of n percolating clusters 
W n H of site and bond percolation on finite square (sq), 
plane triangular (pt), and honeycomb (he) lattices, whose 
aspect ratios approximately have the relative proportions 
1 : \/3/2 : \/3 considered by Langlands, et al. p2| . Using 
nonuniversal metric factors, HLC found that the six per- 
colation models have very nice UFSSF's for E p , P, and 
W n near the critical points H, and at the critical point 
the average number of percolating clusters increases lin- 
early with aspect ratios of the lattices ||. Using Monte 
Carlo simulation, Okabe and Kikuchi found UFSSF's for 
the Binder parameter g jl3j and magnetization distribu- 
tion functions p(m) of the Ising model on planar lattices 
||, and Wang and Hu found UFSSF's for dynamic crit- 
ical phenomena of the Ising model Rj. Based on the 



connection between the g-state bond correlated percola- 
tion model (BCPM) and the g-state Potts model @ Hu, 
Chen, and Lin found UFSSF's for E p and W n of the q- 
state BCPM without using nonuniversal metric factors 
Q. It is of interest to get a deeper understanding of 
the UFSSF's for W n and P for the system with multiple 
percolating clusters. 

In this paper, based on the connection between the 
Ising model, i.e. the two-state Potts model, and the two- 
state BCPM, we use Monte Carlo method to calculate 
the distribution function for the fraction (c) of lattice 
sites in percolating clusters in subgraphs with n perco- 
lating clusters, / n (c), and the distribution function for 
magnetization (m) in subgraphs with n percolating clus- 
ters, p n (m). We find that / n (c) and p n (rn) have very 
good finite-size scaling behavior and they have UFSSF's 
for the model on sq, pt, he lattices when aspect ratios of 
these lattices have the proportions l:\/3/2:\/3. Since W n 
and P of the two-state BCPM for the Ising model may 
be calculated from /«(c), the universality of finite-size 
scaling functions for W n and P are related to universal- 
ity of finite-size scaling functions for f n (c) . The complex 
structure of p(m) for the system with large aspect ratio 
could be understood from the independent orientations 
of two or more percolation clusters in such system. Our 
work suggests many problems for further research. 

The Hamiltonian of the Ising model on an L\ x L 2 lat- 
tice G of Nf, bonds is given by Ti = —JJ2<ij> <J i <J j ~~ 
h Y] ■ Oj, where <7j = ±1, J > and is the ferromag- 
netic coupling constant between the nearest-neighbor 
Ising spins, and h is the external magnetic field. Using 
subgraph expansion, Hu ]T^ | showed that the partition 
function of the Ising model on G may be written as 



Z N = e KN » _ p )^-6(G') 



G'CG 



H [2coBh(£n c (G'))]> 



(1) 



cluster 



where p = 1 - e~ 2K , K = J/(k B T), B = h/{k B T), 
b(G') is the number of occupied bonds in G", and the 
sum is over all subgraphs G' of G, the product extends 
over all clusters in a given G', n c (G') is the number of 
sites in each cluster. When B = 0, Eq. ([!]) reduces to 
the result of Ref. ]l5| . The sites connected by occu- 
pied bonds are in the same cluster; all spins in a clus- 
ter must be in the same direction, which may be up or 
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down. Using Zn of Eq. ([!]), Hu found that the sponta- 
neous magnetization and the magnetic susceptibility of 
the Ising model are related to the percolation probability 
P and the mean cluster size of the (two-state) BCPM 
such that the probability weight for the appearance of a 
subgraph G of b(G') bonds and n(G') clusters is given 
by n(G',p) = p b ( G ">{l -p) Ar ->-fc(G') 2 «(G'). Such connec- 
tion ensures that the phase transition of the Ising model 
is the percolation transition of the BCPM. The exten- 
sion of the probability weight tt(G' ,p) to a g-state Potts 
model is simply achieved by replacing 2 by q in 7r(G",p). 
From n(G',p) we can define 

E p (G, P ) = £ < G »I E (2) 

G' p CG G'CG 

which is called the existence probability of the BCPM. 
Here the sum in the denominator is over all subgraphs 
G' of G and the sum in the numerator is restricted to 
all percolating subgraphs G' p of G. In jl6) , Hu and Chen 
found that Ep(G,p) has very good finite-size scaling be- 
havior. In |9j], Hu, Chen, and Lin found that E p (G,p) 
has UFSSF. 

In the present paper we use the Wolff algorithm |l7]] 
for spin update and study the percolating properties of 
clusters on planar lattices with periodic boundary condi- 
tions for both directions. For the assignment of a bond- 
percolating cluster, we consider free and periodic bound- 
ary conditions in the vertical and horizontal directions, 
respectively; that is, a cluster which extends from the top 
row to the bottom row is a percolating cluster. 

We first consider the fraction of lattice sites in the per- 
colating clusters, c, and denote the probability distribu- 
tion function of c by /(c). The average value of c gives 
the percolation probability P, 

(c) = [ cf(c)dc = P, (3) 
Jo 

and plays a role of order parameter in the percolation 
problem. To study c in subgraphs with exactly n per- 
colating clusters, we decompose /(c) by the number of 
percolating clusters, n; that is, 

oo 

f(c) = J2Uc). (4) 

71=1 

We should note that 

/ f n (c)dc = W n , (n=l,---,oo) (5) 
Jo 

and 

oo 

Y / W n = l-W = E p (G,p). (6) 

n=l 

We may consider the quantity, 



(c) n = / cf n {c)dc, {n = !,■■■, oo), (7) 
Jo 

which is the fraction of lattice sites in the n percolating 
clusters. In other words, we decompose (c) by the number 
of percolating clusters such that 

oo 

71=1 

The probability distribution of the magnetization, 
p(m), is an important quantity in the phase transition 
problem P,|lO|]. Let us decompose p(m) by the number 
of percolating clusters in the same way as in /(c), 

oo 

p(m) = 22p„(m). (9) 

71=0 

Now we have the relation, 




p n (m)dm = W n , (n = 0,- •• ,oo). (10) 



It should be noted that the relation between /„ (c) and 
p n (rn) is not a simple one especially for a system with 
multiple percolating clusters. There are two type of clus- 
ters, that is, the cluster with up spins and that with 
down spins. We may divide the fraction of lattice sites 
in percolating clusters, c, into two classes, c+ and c_. 
By definition, c = c + + c_ and to the leading order, 
m ~ c+ — c_ . If there is only a single percolating cluster, 
m ~ c+ or m ~ c_; thus, m 2 ~ c 2 in the leading con- 
tribution. However, if there are two or more percolating 
clusters, the relation is not trivial, and this is the ori- 
gin of the complex structure of p(jn) for the lattices with 
large aspect ratio a. Therefore, the study of c and m be- 
comes more interesting in the case of multiple percolating 
clusters. 

According to the theory of finite-size scaling if 
a quantity Q has a singularity of the form Q(t) ~ t u 
near the criticality t — 0, then the corresponding quan- 
tity Q(L,t) for the finite system with the linear size 
L has a scaling form Q(L,t) - l v X{tL x l v ), where 
v is the correlation-length exponent and is 1 for two- 
dimensional (2D) Ising model. The finite-size scaling 
is also applicable to the distribution function of Q. 
At the criticality t = 0, we have a finite-size scaling 
form p(Q;L,t = 0) ~ L"^Y(QL U / U ). Thus, we ex- 
pect the following finite-size scaling relations: W n (t) ~ 
X%{tLV% (c) n (t) ~ L-V»Xl{tLV% f n (c;t = 0) ~ 
LP/»Y£(cLP/ u ), p n {m;t = 0) ~ L l3 / u Y^{mL 13 / 1 '), where 
j3 is the order-parameter exponent and is 1/8 for 2D Ising 
model. 

The finite-size scaling functions usually depend on the 
lattice, or other details of the system. However, with 
appropriate choices of nonuniversal metric factors D\ and 
D*, 
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D 2 Q(L,t) = L-"/»X(D 1 tL 1 ^), 
p(Q; L,t = 0) = B 2 L"l v Y(D 2 Qir' v ), 



(11) 
(12) 



the finite-size scaling functions X, Y could become uni- 
versal. This concept of the UFSSF was first proposed by 
Privman and Fisher Q , and has been recently confirmed 
for the percolation problem J5|,|| and for the Ising model 
[§,|]. We should note that the UFSSF's still depend on 
boundary conditions 0. 

To study the finite-size scaling and the universality of 
/„(c), /(c), p n (m), and p{m), we calculate /„(c), /(c), 
Pn(m), p(m), W n , (c) n^ (c), gem and g c for the BCPM on 
sq, pt, and he lattices whose aspect ratios approximately 
have the proportions 1: V3/2: V3 and each kind of lattices 
has two linear dimensions; here g cn — (c)„ 2 /(c 2 )„, g c = 
(c) 2 /(c 2 ), and the second moments of c are defined as in 
Eqs. (0) and ||). We note that g cn and g c have the same 
finite-size scaling property as the Binder parameter p3| . 
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FIG. 1. Universal finite-size scaling functions for the Ising 
clusters on the sq, pt, and he lattices, (a) D^ 1 f n (c)L~ l3 / L ' 
and D^ 1 f(c)L~ >3 ^ l/ at the critical point as a function of 
D 2 cL l3/v , (b) W n as a function of tL 1/v , (c) D 2 {c) n L l3/v and 
D 2 {c)L l3/v as a function of tL 1/v , (d) g cn {= {c) 2 J{c 2 ) n ) and 
9c(= (c) 2 /{c 2 )) as a function of tL 1/v . 

The calculated D^ 1 f n (c)L~ 13 / 1 ' and D^ 1 f(c)L~ 13 / 1 ' at 
the critical point as a function of D 2 cL^l v are shown in 
Fig. |l|(a). The lattice sizes are given within the figure. 
The aspect ratio is a = 4 for sq lattice, and correspond- 
ing equivalent ratios for other lattices. The calculated 
W n , D 2 (c) n LP/» (also D 2 (c)L^ u ), and g cn (also g c ) as 
a function of tL 1 /" (t = (T - T c )/T c ) are presented in 
Fig. ^(b), |l](c), and |l|(d), respectively. The calculated 
D^pJmjL-l 3 /" as a function of D 2 mL^/ 1 ' is shown in 
Fig. ||(a). The metric factors D\ and D 2 for the sq lat- 
tice are chosen as 1 |l8| . The values of Di for the pt and 
he lattices are 1.00 ± 0.01 which are consistent with the 
results of the values of D 2 for the pt and he lattices 
are 1.02 ± 0.01 and 0.98 ± 0.01, respectively. Since we 
have estimated D\ as 1.00 ± 0.01 for the pt and he lat- 
tices, we have omitted D\ in the horizontal axes of the 
figures. Figures |l|(a)-(l|(d) and Fig. [|(a) show that the 
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calculated quantities have very good finite-size scaling 
behavior and the universality is also well satisfied. We 
should note that the metric factors Di and D 2 are the 
same for all quantities. 
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FIG. 2. (a) D 2 ~ 1 p n (m)L- ,3/ ' / at T = T c as a function of 
Di.mLPI 1 ' . (b) p2(m) at T = T c is decomposed into three 
classes, p++(m), p__(m) and (m). 

Figure |^(a) shows that p(m) at T = T c has a broad 
peak centered at m = in addition to two peaks of posi- 
tive and negative to for the system with the aspect ratio 
a = 4 for sq lattice. This is contrast to the case of a = 1 
where p(m) has only two distinct peaks of positive and 
negative to. Such dependence of p(m) on a has already 
been pointed out in Ref. [|oJ. From Fig. §(a), we see 
that the broad peak of p(m) centered at m = mainly 
comes from p 2 {jn). There are two types of Ising clus- 
ters, that is, the clusters with up spins or the clusters 
with down spins. Therefore, if there are many percolat- 
ing clusters, the combination of the percolating clusters 
with up spins and those with down spins makes it possi- 
ble that the total magnetization becomes close to 0. It is 
known that the normalized fourth moment of to, or the 



Binder parameter, at the critical point depends on the 
aspect ratio (l(]jnj • The origin of such a dependence can 
be attributed to the structure of many percolating clus- 
ters. To clarify this situation, we decompose pziyn) at 

T = T c into three classes, p ++ (m), p__(m) and p_| (m) 

shown in Fig. ^(b). We assign three peaks in p 2 (to) by the 
contribution from p ++ (to), p (to) and p_| (to). Exam- 
ples of snapshots of the Ising system with two percolating 
clusters are presented in Figs. ||(a) and||(b). In Fig. S(a) 
both percolating clusters are up; in Fig. ||(b) one per- 
colation cluster is up and another percolating cluster is 
down. 
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FIG. 3. Examples of snapshots of the Ising system with 
the aspect ratio a — 4. Up and down spins are represented 
by open and closed circles, and active bonds are represented 
by solid line. Percolating clusters are distinguished by shaded 
area. There are two percolating clusters with up spins in (a), 
whereas one percolating cluster is up and the other is down 
in (b). 

From W n , we may calculate the average number of 
percolating cluster by (n) = J2 n nW n . At the critical 
point, the values of W n and (n) as a function of the as- 
pect ratio a = Li/L 2 are plotted in Fig. [|(a) and ||(b), 
respectively. We see that (n) increases linearly with a 
for large a, which is similar to the case of random per- 
colation H^O). The slope of (n) versus a in Fig. ||(b) is 
approximately 0.5. 

Following the study of W n for random percolation by 
Hu and Lin ||, there have been many analytic and nu- 
merical studies of W n in different random percolation 
problems pi] ]; it is of interest to extend such studies 
to the BCPM of the Ising model. On the other hand, 
we can extend the study f n (c) and (c) n to the random 
percolation problem. It is interesting to compare the 
results for the BCPM and those for the random perco- 
lation problem. We may also extend the present study 
to the bond-diluted or the site-diluted Ising model which 
can be mapped into percolation models |22(] . The criti- 
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cal phenomena of the percolating properties of the Ising 
model are governed by the Ising fixed point (for example, 
v = 1), whereas at the percolation threshold the critical 
phenomena are governed by the random percolation fixed 
point (y = 4/3). The crossover from the Ising fixed point 
to the random percolation fixed point in the process of 
dilution is highly interesting, especially for the proper- 
ties depending on the number of the percolating clusters. 
The studies in these directions are in progress. 





FIG. 4. (a) W n at T = T c as a function of a = Li/L 2 . (b) 
(n) at T = T c as a function of a. 
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